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Abstract 

We study a simple dynamical model exhibiting sequential dynamics. We show that 
in this model there exist sets of parameter values for which a cyclic chain of saddle 
equilibria, Ofc, k = have two dimensional unstable manifolds that contain 

orbits connecting each Ok to the next two equilibrium points Ok+i and Ok +2 in the chain 
(Op+i = Oi). We show that the union of these equilibria and their unstable manifolds 
form a 2-dimensional surface with boundary that is homeomorphic to a cylinder if p is 
even and a Mobius strip if p is odd. If, further, each equilibrium in the chain satisfies 
a condition called “dissipativity,” then this surface is asymptotically stable. 


1 Background 

In the last decade it became clear that typical processes in many neural and cognitive 
networks are realized in the form of sequential dynamics (see m, m, m, m, m, 
[2] and references therein). The dynamics are exhibited as sequential switching among 
metastable states, each of which represents a collection of simultaneously activated nodes 
in the network, so that at most instants of time a single state is activated. Such dynamics 
are consistent with the winner-less competition principle In the phase space of 

a mathematical model of such a system each state corresponds to an invariant saddle 
set and the switchings are determined by trajectories joining these invariant sets. In the 
simplest case these invariant sets may be saddle equilibrium points coupled by heteroclinic 
trajectories, and they form a heteroclinic sequence (HS). This sequence can be stable if all 
saddle equilibria have one-dimensional unstable manifolds [2], in the sense that there is an 
open set of initial points such that trajectories going through them follow the heteroclinic 
ones in the HS, or unstable if some of the unstable manifolds are two-or-more dimensional. 
In the latter case properties of trajectories in a neighborhood of the HS were studied in 
[2] and [3]. General results on the stability of heteroclinic sequences have been obtained 
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by Krupa and Melbourne [HI US] in the form of necessary conditions that may also be 
sufficient in the presence of certain algebraic conditions. 

An instability of a HS may be caused if some initial conditions follow trajectories on the 
unstable manifold different from the heteroclinic ones. M. Rabinovich suggested [22] con¬ 
sidering the case when all trajetories on the unstable manifold of any saddle in a HS are 
heteroclinic to saddles in the HS, which implies an assumption that the HS is, in fact, a 
heteroclinic cycle. In this case one can expect some kind of stability, not of the HS, of 
course, but of the object in the phase space formed by all heteroclinic trajectories of the 
saddles in the HS. We study in our paper the Rabinovich problem. We deal here with 
the generalized Lotka-Volterra model [2] that is a basic model of sequential dynamics for 
which unstable sets are realized as saddle equilibrium points. All variables and parameters 
may take only nonnegative values, so we work in the positive orthant of the phase space 
M”. We impose some restrictions on parameters under which all unstable manifolds of the 
saddle point are two-dimensional and all trajectories on them (in the positive orthant) are 
heteroclinic in some specific way (see below). We prove that they form a piece-wise smooth 
manifold homeomorphic to the cylinder if the number of the saddle points in the HS is even 
or to the Mobius band if it is odd. We prove also that under the additional assumption 
that each equilibrium is dissipative (see below), then this manifold is the maximal attractor 
for some absorbing region (in the positive orthant). Trajectories in this region may follow 
different heteroclinic trajectories and may manifest some kind of weak complex behavior. 

Although our motivations are neurological, we observe that heteroclinic networks (and 
thus, potentially, high-dimensional counterparts of the same) are ubiquitous, appearing in 
applications that range from celestial dynamics m to evolutionary game theory P. 


2 Notation and Results 


We begin our consideration of two-dimensional heteroclinic channels with the study of a 
series of Lotka-Volterra equations. Lotka-Volterra models are widely used in the context of 
heteroclinic sequences where all the unstable manifolds are one dimensional, so that each 
equilibrium is connected to exactly one subsequent equilibrium (e.g. [laiiHiiinii). It has 
been recently shown that they ask provide a general method for embedding directed graphs 
into a system of ordinary differential equations [ 6 |. We remark that although [ 6 ] allows 
graphs of high valency to be modeled by heteroclinic networks, and some work has been 
done on the stability of such systems (e.g. US], which considers the competing dynamics 
of the “overlapping” channels Ai —>• A 2 —>• A 3 —>• Ai and Ai —)• A 2 —)• A 4 —)• Ai), study 

of heteroclinic networks has usually only considered one-dimensional unstable manifolds. 
In keeping with the discussion of the introduction, we consider the dynamics of a system 
such that initial conditions arbitrarily close to any of p < n saddle nodes may be mapped 
into neighborhoods of either one of two other saddle nodes. In particular, a simple general 
model for heteroclinic sequential dynamics was given in ID by 

n 

Xi = Fi{x) = Xi{ai-'^PijXj) for i = 1, ...,n, (1) 

i=i 
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Figure 1: A rough representation of the dynamics of the system defined by Equation [H 
with p = 5. 


where all of the parameters are assumed to be positive. The constants pij have biological 
meaning, representing inhibition of mode i by mode j. For the sake of simplicity, we assume 
that Pa = 1 for all i. Further, since the variable x is assumed to encode biological informa¬ 
tion that is necessarily non-negative, e.g. activation levels or chemical concentrations, we 
restrict the system to the hrst closed orthant, = {x G M” : x* > 0 , 1 < z < n}. The 
system is so constructed that it contains saddle points lying on the axes, with ai being the 
z-th coordinate of the z-th saddle along the z-th axis > 0), i.e., the system ([1]) has n 
equilibrium points of the form Ok = ( 0 , ..., 0 , ( 7 ^, 0 , ...), for k = 1, ..,n. There may be other 
equilibria as well, but they are not relevant for our purposes; we only study transitions 
between the n equilibria just dehned. 


In the present work we suppose that the first p equilibria points are sequentially connected 
by a set of 2-dimensional unstable manifolds. For each k, 1 < k < p, there will be a 
heteroclinic orbit connecting Ok to Ofc+i and a heteroclinic orbit connecting Ok to Ofc+ 2 - 
Furthermore, the system is closed in the sense that Op+i = Oj, i.e. p is the modulus of 
the subscript. In the following, we will consider the restrictions necessary to enforce such 
dynamics. 


To ensure that there are heteroclinic trajectories between Ok and both Ok+i and 0^+2) 
we apply eigenvalue conditions to the system. Consider first Oi. The linearization of the 
vector held F(x) (HJ at Oi is given by the upper triangular matrix: 



/ 

-£7l 

— 0'lPl2 

—<7i/513 

— <7i/?14 



0 

<72 — P2lO'l 

0 

0 

DF{Oi) = 


0 

0 

0 

0 

<73 — /531<7i 

0 

0 


<74 — P41<Ti 






0 


1 

0 

0 

0 
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^iPln 

0 

0 


( 2 ) 


Pnl^l j 


and so the eigenvalues appear on the diagonal. The matrices DF{Ok) have a similar simple 
structure, zeros everywhere except on the diagonal and on the k-ih. row. It is then easy to 
see that the eigenvalues of DF at Ok are 

and = aj - PjkfJk, for j / k. 

In particular, the eigenvalues are all real. One can also see from the structure of ([2]) that 
DF{Ok) has a full set of eigenvectors, even if some eigenvalues are repeated. 
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Note that because of the particular form of the equations, all coordinate axes, planes and 
hyperplanes are invariant. Thus, in order that trajectories can travel from Ok to Ofc+i or 
, it IS sufficient to put the restriction 0 ^ Afc+ 2 >^fc+i> and < 0 otherwise, to ensure 

that they can only go in those directions. In particular, for each k, 1 < k < p, we require 
that 


0 < min{cTfc+i - pk+i,kc^k} (indices mod p), (3) 

2 = 1,2 

~ ~ Pjk^k <0, for j 7 ^ /c. A; + 1, A: + 2 mod p. (4) 

Note also that —ak < 0. These inequalities guarantee that each equilibrium is a hyperbolic 
saddle with 2 unstable directions and n — 2 stable directions. Let = W‘^{Ok) n M” be 
the unstable manifold of Ok restricted to the positive orthant. We show below that 

T^UiW^UOk) 

k=l 

forms a piecewise smooth surface that we will classify topologically as follows. 

Theorem 2.1. Suppose that inequalities Q and ([3]) hold for each k, 1 < k < p and that 
each unstable manifold fT“(Ofc) is contained in a compact forward invariant set as specified 
in Lemma \S.10\ (see also Remark 3.11). When p is even, the union of unstable manifolds 
r is homeomorphic to a cylinder. When p is odd, T is homeomorphic to a Mobius strip 


The proof of this theorem, which is slightly involved, is left for the appendix. 

Consider the following definition. 

Definition 2.2. Let and be the set of stable and unstable eigenvalues, respectively, 
of the linearization of a vector field at a saddle equilibrium, i.e., maxiie(S^) < 0 and 
mini2e(S“) > 0. We say that the saddle is dissipative if 

maxi?e(S“) < — maxi?e(S^). 

In other words, the weakest stable eigenvalue is stronger than the strongest unstable eigen¬ 
value. (See w 

In terms of the specific vector field under study all the eigenvalues in question are real and 
for each k we have: 


max{(Tfc+j - < min - p,fc(Tfc|, o-fc} (indices mod p). (5) 

i=l,2 jjik,k+l,k+2 


The main goal of this manuscript is to show that, under the condition that each saddle 
equilibrium is dissipative, then T is asymptotically stable. Specifically, our main theorem 
is: 

Theorem 2.3. Suppose that inequalities (ED, ED and ED hold for each k, 1 < k < p and 
that each unstable manifold W'^{Ok) is contained in a compact forward invariant set as 
specified in Section\3f3 Then T is asymptotically stable. 
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The proof of Theorem 12.31 breaks roughly into two independent pieces. We start by consid¬ 
ering the trajectory of a representative point that is e-close to T, but is distant from each 
of the fixed points Oj. In such a case, the dynamics of the system are controlled largely 
by three consecutive saddles Oj, Oi+i, and Oj+ 2 . In Section [3l we consider the restriction 
of ([I]) to three consecutive dimensions; we will gain information on the full-dimensional 
system by viewing it as a perturbation of this restriction. The second part of the proof 
is to consider the dynamics as a trajectory passes near a hxed point; we consider this in 
Section [H 

In Section 02] we prove the main theorem. In Section [4.31 we show that there is a non-empty 
parameter set for which the conditions of the theorem are satisfied. 


3 Three Dimensions 


We begin our proof of Theorems 2.1 and 2.3 with a study of the restriction of the system 
to three dimensional sub-spaces corresponding to three consecutive coordinate directions. 
Through a series of geometric lemmas, we prove the main result of the section. Theo¬ 
rem 13.121 which provides information on the behavior of trajectories inside this invariant 
subspace. This will be used in later sections to complete the proofs of the main results by 
providing information on those parts of the full phase space where all but three coordinates 
are small. 


3.1 Set-up 

Let Oj, Oj+i, and Oi +2 be any three consecutive equilibria and restrict the system ([T|) to the 
three dimensions spanned by consecutive coordinates Xi, Xj+i, and Xj+ 2 . For convenience, 
we will refer to these three variables as xi, X 2 , and X 3 . Our goal in three dimensions is to 
show the existence of a compact, forward invariant set containing Oi, O 2 , and O 3 such that 
any trajectory with an initial value in the interior of that set converges to O 3 . 

Restricted to three dimensions, the equations ([I|) are reduced to 

xi = xi{ai -xi- P 12 X 2 - P 13 X 3 ), 

X 2 = X 2 {(J 2 -X 2 -P 21 X 1 -P 23 X 3 ), (6) 

X 3 = X3{a3 -X 3 - p 3 iXi - P 32 X 2 ). 

The restriction of the dimension of the unstable manifolds via the eigenvalue conditions, 
and the positivity conditions on (Ti,(T 2 ,and <73 yield the following inequalities: 


fJi < 0 < Uj - Pjiai, j = 2, 3, 

(7) 

—172 < 0 < (73 — P320-2, 

( 8 ) 

<71 — P12<72 < 0, 

(9) 

<7l — Pl3<73 < 0, 

( 10 ) 

<72 — /023<73 < 0. 

( 11 ) 
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Of those inequalities, ([7|) controls the behavior of the system at Oi, ([ 8 ]) - Q control the 
behavior of the system at O 2 , and (fTOl) - (fTT]l at O 3 . The point ((Ti, 0,0) is a saddle with 
a two-dimensional unstable manifold, ( 0 ,( 72 , 0 ) is a saddle with a one-dimensional unstable 
manifold, and (0, 0 ,( 73 ) is a sink for the system Q. We remark that although ([ 6 |) has the 
form of the May-Leonard model, the particular parameter restrictions under consideration 
yield simple dynamics (see Theorem I3.12p , and prevent the more complex behavior usually 
studied in that context. In particular, they are inconsistent with the symmetric May- 
Leonard model as it was introduced in [El- 



Figure 2: Illustration of dynamics of the system projected onto the first three coordinates. 


All consecutive triplets (xj,Xj+i,Xj+ 2 ) (where Xp+j = Xj) of coordinates possess the same 
qualitative dynamics. 

For each z, we will be interested in the points where Xj = 0; each such set is the union of 
the plane Xj = 0 and some “nontrivial” plane. We designate those planes: 


Pi ■= Wi - P12X2 - P13XS - xi = 0 }, 

P2 ■= {cr2 - P21X1 - P23X3 - X2 = 0}, 

P3 ■= {0-3 - P3lXl - P 32 X 2 -X3 = 0}, 


( 12 ) 

(13) 

(14) 


where Xj = 0 on Pj. 

We will also refer to the plane passing through the points ((Ti, 0,0), (0, ( 72 ,0), and (0,0, ( 73 ), 
which we denote by S. Observe that S is given by the equation 



CJi (72 (73 


For ease of discussion, we will also name the coordinate planes: 


P12 ■= {(xi,X2,X3) : X3 = 0}, 
P23 ■= {{xi,X2,X3) : xi = 0}, 
-P 13 := {(3:i,3:2,3:3) : X 2 = 0 }. 
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Because of the positivity conditions on the parameters and variables, we may use Pi, P 2 , 
P 3 , and S as shorthand for the intersection of those planes with the first octant without 
the risk of confusion. We observe that the intersection of S and of each Pi with is a 
compact triangle, a fact we will use repeatedly in the following section. 

In order to describe the dynamics of orbits, we are interested in when one plane lies “above” 
another in the first octant. Consider the following definition: 


Definition 3.1. Observe that each plane Pi can be written as the graph of a function 
P(xi, X 2 ). A plane Pi dominates a plane Pj if {xi,X 2 , P(xi, X 2 )) G Pi and (xi, X 2 , P(xi, X 2 )) G 
Pj implies that z^{xi,X 2 ) < z^{xi,X 2 ) for all xi,X 2 > 0. d plane Pi is dominated by a 
plane Pj if Pj dominates Pi . 


3.2 Geometrical Lemmas 

Each plane Pi divides into two regions, one where Xi is positive and another where it is 
negative. This allows information about Xi to be gained from purely geometric information. 
For example, xi is positive below Pi, and negative above it. Since P 3 dominates Pi (i.e. 
is always above it), we instantly see that xi\p^ < 0 (see Corollary 13.51 given below). We 
introduce geometric lemmas giving the information we can gain in this way; the proofs of 
all of them are parallel to one another, and can be summarized as follows: since the system 
is restricted to the first octant, we compare two planes (triangles) by seeing where they 
intersect the xi, X 2 , and X 3 axes. Denote the compact triangle thus formed by a plane S 
as Ts and the non-zero component of its vertex on the i-th axis as 5*. Then a plane S 
dominates a plane P if P* < S'® for i = 1,2,3, with at least one of those a strict equality, 
i.e. if its vertices are farther from the origin. 

Lemma 3.2. The plane Pi is dominated by the plane S. 


Proof. We consider where each plane intersects each axis: 


• The planes Pi and S both intersect the xi—axis at the point Oi. 

• The plane Pi intersects the X 2 axis at (0,^,0), while T, intersects the axis at O 2 . 
We know from Q that < <^ 2 - 

• The plane Pi intersects the X 3 — axis at (0,0, ^), while S intersects the axis at O 3 . 
We know from (fTOjl . that ^ < <73- 

Since P( < S® for all i. Pi is dominated by S. □ 

Lemma 3.3. The plane S is dominated by P 3 . 


Proof. We consider where each plane intersects each axis: 
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• The plane P 3 intersects the xi—axis at the point 0; 0); while S intersects the axis 

CTq 

at Oi- We know from ([7]) that cJi < —. 

/531 

• The plane P 3 intersects the X 2 axis at (0,^,0), while S intersects the axis at O 2 . 
We know from (l 8 |) that ^'2 < 

• The planes P 3 and S both intersect the X 3 axis at O 3 . 

Since S* < P 3 for all i, P 3 dominates S. □ 


The property “is dominated by” is clearly transitive, so the following corollary holds. 
Corollary 3.4. Pi is dominated by P 3 . 

Since xi < 0 above the Pi plane, the following corollary follows immediately. 

Corollary 3.5. On the plane P 3 , xi < 0. 

If we further had that S dominates P 2 , then the eigenvalue conditions introduced in [4] and 
summarized as (3) - ( 8 ) would be sufficient to ensure the existence of a positively invariant 
region. It happens, however, that this is not the case. 

Lemma 3.6. The plane S neither dominates nor is dominated by P 2 . 

Proof. We consider where each plane intersects each axis: 

• P 2 intersects the xi—axis at the point (^, 0 , 0 ), while S intersects the axis at Oi. 

We know from ([7]) that and therefore S does not dominate P 2 . 

• P 2 intersects the X 3 — axis at (0,0, ^), while S intersects the axis at O 3 . We know 
from (fTT]l that ^ < <^ 3 ) and therefore P 2 does not dominate S. 

Thus, neither plane dominates the other. □ 


The situation is somewhat salvaged by the following. 

Lemma 3.7. If — < —, then P 2 is dominated by P3. 

Proof. In the proof of Lemma 13.31 (second bullet point), we established that P| < P|. 
In the proof of Lemma 13.61 (second bullet point), we established that Pf < P 3 . All that 
remains for P 3 to dominate Po is for P^^ < Pn, which occurs if and only if — < □ 

The parameter restriction in (j3.7jl written in terms of the general systems gives that for 
each fc, 1 < < p, 

o'fc+i ^ o'fc +2 (indices mod p). (15) 

Pk-\-l^k Pk+2,k 

Similarly to Corollary 13.51 we have the following: 
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Corollary 3.8. In the region of parameter space where the hypotheses of Lemma [X 
satisfied, X 2 \p,^ < 0 . 


are 


We note here one additional observation. 

Lemma 3.9. The rectangular box: 

B = {x : 0 < Xi < ai,i = 1,... ,n} 

is forward invariant with respect to the system (HI)- Further, the unstable manifolds W]f are 
all contained in B. 


Forward invariance follows immediately from the differential equations ([T]) and the assump¬ 
tion that pii = 1. The conclusion that the unstable manifold at Oi is inside B follows easily 
by noting that the unstable eigenspace at Oi (restricted to the first orthant) is strictly 
inside B. 



Figure 3; The plane Pi (dotted) is dominated by S (solid) which is dominated by P 3 
(dashes). If the hypotheses of Lemma [3^1 are satisfied, then the plane P 2 (dashes and dots) 
is dominated by P 3 . 

Rather than requiring that P 3 dominates P 2 , we could require only that P 3 dominates P 2 
inside the forward invariant box B. Since P 3 is strictly outside of B by Lemma 13.31 (see 
Figure [ 3 ^ . this produces a strictly greater set of allowable parameter values. 


3.3 Existence of a Positively Invariant Region 


We remarked that our goal in three dimensions was to show the existence of a compact, 
forward-invariant set whose trajectories converge to O 3 . We now carry this out. 


Lemma 3.10. In the region of parameter space where the inequalities (0) - m and the 
hypotheses of Lemma \3. 7| are satisfied, the planes P 3 , P 23 , P 13 , and P 12 enclose a positively 
invariant region, i.e. no trajectory leaves in positive time. 
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Proof. No trajectory can leave through any of the Xj = 0 planes, since Xilx^-piane = 0 . The 
outward normal vector to P3 is N = (/O31,/932, 1 ), and its scalar product with the vector 
field ([H) is given by h • F = psixi + P32X2 when X3 = 0 . This is negative by Corollary 13.51 
and Corollary 13.81 □ 


Remark 3.11. Since the pair of inequalities psixi < 0 and P32X2 < 0 form a suffieient, 
but not neeessary, condition for the inequality p3ixi + P32X2 < 0 to be satisfied, a positively 
invariant region might exist even if the hypothesis of Lemma 3.1 is not satisfied. For 
instance, since the box B (Lemma \ 3 .ff\) is forward invariant in may be that X2|p3ns is 
negative even when the hypotheses of Lemma \ 3 . 7 \ fail. Lemma [ 3 . 10 \ can therefore be extended 
to a region of parameter space that includes the region defined by Lemma \ 3 . 7 | as a proper 
subset. 


Theorem 3.12. Suppose that P 3 , xi = 0 , X2 = 0, and X3 = 0 enclose a positively invariant 
region. Any trajectory (ff{x) in the above-described region that is not contained in P12 (the 
3^3 = 0 plane) goes to (0,0,0-3) as t ^ +00. 


Proof. Any trajectory in a compact positively invariant region has a non-empty cj-limit 
set. Fix a trajectory with initial condition x in the interior of the region, and let q be an 
arbitrary point in the w-limit set Llx- The proof breaks into three parts: we prove that q 
lies on P3, that it lies on P23) and that it lies on P13; the intersection of these planes is 
( 0 , 0 ,(T3). The proofs of the second and third statements are essentially the same as the 
proof of the first statement, which we treat in detail. 

By way of contradiction, suppose that q does not lie on P3. First of all, note that q does 
not lie on P12, because by assumption, the x ^ P12, and since the trajectory increases in 
the X3 variable, (j)^{x) cannot approach P12 = {xs = 0} 

Since X3 is continuous, and q lies on neither P3 nor P12, the regions where X3 = 0 , we can 
find a spherical neighborhood Or{q) centered at q with radius r such that X310^.(9) > e > 0 
for some e. It follows that q is not a fixed point, and Llx is not the singleton {g}, since 
w—limit sets are forward invariant |21] . 

The X3 component of ((f{x) is non-decreasing in the positively invariant set and it is bounded 
above, since it is bounded above by <73. Thus it has a limit X3. It follows by continuity that 
the X3 component of any point in Fix is also Xg. We have supposed that q G Fix but q ^ P3] 
now consider (p^{q). Since X3 > 0 at q, and thus in a neighborhood of q, the X3 component 
of (f*^{q) must be strictly increasing as a function of time along the forward solution near q. 
This contradicts that the X 3 component is Xg everwhere on Fix. 

We repeat the argument twice. First, q must lie on the plane P23; otherwise xi{q) < 0, 
and the same contradiction will be obtained. Then, using the same argument, we see that 
it must lie on Pl 3 - We observe here, with reference to Lemma 13.71 that even if P2 is not 
dominated by P3, it is dominated by it on the restriction to P23) and the argument therefore 
goes through. Thus q must lie on the single point, (0,0, fig), and our proof is complete. □ 
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The Dynamics in the Full Phase Space 


4.1 Local Dynamics Near an Eqnilibrium Ok 

We assumed in ([5l) that the saddle points Ok are dissipative. Denote by v the ratio: 

min|{i?e(S")}| 

max{i2e(S“)} 

We call V the minimal saddle value of the equilibrium (see [2Z])- The equilibrium is dissi¬ 
pative if 1 / > 1. 

The main implication of this assumption is that a trajectory starting at an initial value 
at a distance e from the stable manifold of the saddle comes to a point of distance on the 
order of e'^ from the unstable manifold after going through a neighborhood of the saddle. 
Formulating this result more strictly, we label the variables in a neighborhood of Ok into 
the 2-dimensional unstable subspace ry = {xk+i,Xk+ 2 ) and the (n — 2)-dimensional stable 
subspace ^ = (xj),j ^ k+1, k+2. Let | • | denote the sup norm in these local coordinates. By 
the Stable Manifold Theorem, for each k there exists 6 k > 0 such that in a (5fc-neighborhood 
of Ok, the unstable manifold W^{Ok) is the graph of a (smooth) function ^ = h'^{rj). Let 6 
be the minimum of these 6 k and consider the 5-neighborhood 14 of each Ok, k = 1,... ,p 

For fixed k, and 0 < sufficiently small, define a pair of sections, 

'S'o = {(C,?/) : l?l = hi < efc} and 5i = {(C, ?y) : |^| < 5, |ry| = 5}. 

By the classical Shil’nikov variables technique (see m), there exists €k sufficiently small 
so that every forward solution starting at Sq will intersect before leaving 14. Let e be 
the minimum of the e^’s needed in the neighborhood of each Ok- Let T denote the time at 
which such a solution intersects Si. 

Theorem 4.1. If 6 and e are sufficiently small and if {f,{0),ri{0)) G Sq and {f,{T),r]{T)) G 
Si then h(T)| < C\r]{0)\'^~^, where C > 0 is independent of the initial point and n — e > 1. 


The sections Sq and Si , together with Theorem 14.11 are illustrated in two dimensions in 
Figure im 


Proof. For ease of notation, consider Oi- Note from ([2]) that the eigenvector corresponding 
to the first diagonal element, —ui, can have only one non-zero component and that is in 
the xi direction. The eigenvector corresponding to the j-th diagonal element, j 1, has 
non-zero components in the xi and Xj directions only. 

Now take into acconnt that the hyperplane {(x 2 = X 3 = 0)} is invariant under the flow of 
the equations and tangent to the stable eigenspace of DF{Oi). It thus coincides locally 
with and contains W^{Oi). 

Next note that the unstable eigenspace has no non-zero components in the coordinate 
directions X 4 ,X 5 ,... ,Xn and that the hyperplane where they are zero is invariant under 
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Figure 4 : A schematic diagram of Theorem 14.11 Initial conditions on a section intersecting 
a stable manifold are mapped in finite time to a section intersecting an unstable manifold, 
resulting in a contraction. 

the flow. It thus follows that the X4,X5,... ,Xn coordinates of the unstable manifold are 
all zero. Thus the unstable manifold W^{Oi) is given locally in the d-neighborhood of Oi 
as the graph xi = ai + /i“(x2,X3) where /i( 0 , 0 ) = 0 and h is smooth. Consider the local 
change of variables that “straightens out” the unstable manifold: 

Xi = xi - (Ji - h{x2,x^). ( 17 ) 

Under this smooth change of coordinates ([T]) becomes: 

Xi = -aiXi - /(Ai, X2, ...,Xn) 

±2 = X2 {(72 -X2- P2l{Xl + CTl + h{x2,Xz)) - P 23 X 3 “ ^243^4- P2nXn) 

X3 = X3 (0-3 - X3- P3l{Xi + CJi + h{x2,X3)) - P 32 X 2 - P 34 :X 4 - PSnXn) 

X4 = X4 (0-4 -X4- P4liXi + (Ti + h{x2,X3)) - P 42 X 2 - P43X3 - PAnXn) 


Xn — Xn {(7n Xn Pnl{X 4 + (7i + h{x 2 , X 3 )) Pn 2 X 2 Pn3X3 ' ' ' Pn,n—lXn—l) ■ 

Since ly^(Oi) is contained in the coordinate hyperplane, we do not need a change of vari¬ 
ables corresponding to W^. Noting the /( 0 ,X 2 ,... ,a;„) = 0, we use the MVT to define a 
new function _ j^^Xi,X 2 , ■ ■ ■, Xn), and rewrite the equations to obtain: 

Xi = —aiXi + fi{Xi,X 2 , ..., Xn)Xi 

±2 = X 2 {(72 - P 2 l( 7 l -X 2 - P 2 l{Xl + h{x 2 ,X 3 )) - P 23 X 3 - P 24 X 4 - P2nXn) 

X 3 = X 3 ( 0-3 - P 3 iai -X 3 - P 3 l{Xi + h{x 2 ,X 3 )) - P32X2 - P 3 AX 4 - P3nXn) (19) 

3:4 = 3:4 (<74 — P4iai — X4 — P4l{Xi + h{x2,X3)) — P42X2 — P43X3 — ■ ■ ■ — P 4 nXn) 
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Now note that distances from W^{Oi) and W'^(Oi) are not effected by the coordinate 
change (fT71) . If we now denote ^ = (Xi, X 4 ,..., x„) then |,^| is the distance to the unstable 
manifold and \r]\ is the distance to the stable manifold (in the sup norm). 


It now follows immediately from (|19p that inside the 5 neighborhood of Oi the unstable 
directions satisfy the estimates: 

X 2 < X 2 (o -2 -/ 02 i 0 'i - e“) and X 3 < X 3 (cj 3 - psicxi - e“), 

where we can take e“ > 0 arbitrarily small by starting with 6 and e small. Thus, by a 
simple application of Gronwall’s inequality, solutions starting on Sq and remaining in our 
6 neighborhood of Oi must satisfy: 

where is the maximal unstable eigenvalue, i.e.. 




- pjiai}. 
3 — 


Thus if T is defined by \r]{T)\ = 6, then 


T > 


1 

Xuu _ 



Similarly, using Gronwall’s inequality once again we obtain: 

\m\ < 

where A^^ is the “leading” stable eigenvalue, i.e. since the eigenvalues are real, the negative 
eigenvalue with the smallest absolute value. In terms of our parameters: 

A^* = max{-fji, {aj - PjiaiJ = 4,... ,n}}. 


Thus 

|{(r)| < c(«)|„(0)|(>'-'' 

where e can be taken arbitrarily small. By the assumption that Oi is dissipative, v > 1 
and so we can make v — e > 1. 


□ 


4.2 Stable Heteroclinic Surface 

We are now in a position to prove that the union of the unstable manifolds of our system, 
restricted to the positive orthant, form an asymptotically stable forward invariant set under 
appropriate parameter restrictions. 

Theorem fTheorem l2.3l restated). Suppose that inequalities ([3]), dH) and ([5]) hold for each 
k, 1 < k < p and that each unstable manifold W'^{Ok) is contained in a compact forward 
invariant set as in Section Then T = ULi(>r? U Ok) is asymptotically stable. 
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We consider the role each inequality of the hypothesis plays in the theorem. The inequalities 
([3]) and dH) state that 

0 < min{(Tfc_|_j — Pk+i,k(^k\ (indices mod p), for each 1 < A; < p, and 
2 = 1,2 ’ 

~ ~ Pjk^k <0, for j 7 ^ A:, A: + 1, A: + 2 mod p, for each 1 < A: < p. 


These inequalities are fundamental to the problem we are considering. They ensure that 
there is a heteroclinic channel from each equilibrium Ok to the equilibria Ok+i and Ok +2 
(the hrst inequality), and that there are no heteroclinic channels to the other equilibria 
(the second inequality). 

The inequality ([5]) states that for each k, 

max{ak+i - Pk+i,kcrk} < min {\(7j - pjkcrk\,ak} (indices mod p). 

2 =1,2 j^k,k^l,k-\-2 

This is the dissipativity condition. Informally, it may be taken to say that when a trajectory 
that is close to one of the unstable manifolds whose union is T passes near one of the saddle 
fixed points, the distance of the trajectory to T contracts exponentially. 

It will become clear, from the proof of Theorem 12.31 that ([5]) is a much stronger condi¬ 
tion than is necessary. It ensures not only asymptotic stability, but a sort of monotonic 
asymptotic stability, such that whenever the trajectory passes near some Ok, its distance 
to r contracts exponentially. If, on the other hand, Q held for some, but not all, values 
of k, then the trajectory would at times contract exponentially towards T, and at other 
times drift away from T, and stability would depend on how these attractive and repulsive 
forces average over time. Attempting to formulate a replacement condition for ([5]) that is 
necessary as well as sufficient is extremely nontrivial. 

The final hypothesis, that each unstable manifold is contained in a compact forward in¬ 
variant set, is necessary. We have seen one set of inequalities that ensure that such sets 
exist, 

(indices mod p), 

Pk-\-l,k P/c+2,fc 

which are sufficient but not necessary. In Section [4.31 we see that there is a nonempty open 
region of parameter space where the hypotheses of Theorem 12.31 hold. From our current 
discussion, we see that the theorem applies to a larger region. 

We recall that T is asymptotically stable if given any neighborhood 17 of T (restricted to 
M”) there exists an e-neighborhood of T, say ^(r) C K”, such that if xq G I7;(r) then 
x{t,xo) G [/ for t > 0 and limt^oo dist(a;(t, xq), T) = 0 , where x(t, xq) is the solution of the 
initial value problem ([1]) with x(0,xo) = xq. When we speak of an open e-neighborhood of 
r, we are speaking of a set that is open in the subspace topology. That is, an e-neighborhood 
in M” is an e-neighbor hood in M” intersected with M”. 


Proof of Theorem \2.3[ For each k = l,...,p, let V{Ok) be a sufficiently small 6 neighbor¬ 
hood of Ok, such that Theorem [JT] can be applied within each V{Ok) and 6 does not depend 
on k. 
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Let zq be a representative point at an initial condition e-close to L. Choose e such that 
€< 5. Then we can classify the dynamics as either local if zq G V{Ok) for some k, and 
global otherwise. 

We observe that if zq = Ok for any k, its behavior is trivial, and likewise, if zq lies on a 
coordinate plane, it remains on that plane while converging exponentially to some Ok- We 
therefore assume without loss of generality that neither of these cases hold. 

Suppose that zq ^ V{Ok) for any k. The point is e-close to T, and since T is a hnite 
union, we can say that zq is e-close to where a is fixed and depends on zq. Consider 
the projection of the system onto the three-dimensional subspace spanned by the axes 
Xa+i, and Xa+ 2 - In three dimensions, the specific route a solution takes has not been 
important; a trajectory in the invariant set may go straight to a neighborhood of 0 ^+ 2 j or 
it may detour to Oa+i, but the net result is the same iTheorem 13.121) . We now formally 
differentiate between these two cases. 

We consider two cases: either the positive semitrajectory of zq intersects V(0a+2) without 
first intersecting V{Oa+i) (case (i)), or the positive semitrajectory intersects V{Oa+i), then 
intersects V(Oa+ 2 ) (case (ii)). 

Before proceeding, we recall the definitions of Sq and Si given in Theorem l4.ll and similarly 
define such sections Sq, 5^ for 1 < < p. Without loss of generality, we assume that 2o € •S'l. 

Suppose that case (i) occurs. For each k = 1, ...,p, let V{Ok) be the projection of V{Ok) 
into the three dimensions spanned by Xk-,Xk+i-, and Xk+ 2 '-, note that Xk is negligable for 
/c 7 ^ a, a-|-1, and a + 2. Then the projection of the orbit onto intersects V{Oa+ 2 ) before 
it can intersect V{Oa+i)- In we know that all trajectories inside of that do not 
intersect V (Fh+i) come to a neighborhood of 0^+2 in bounded time, where the bound does 
not depend on the initial condition. We may consider the non-projected, full-dimensional 
space as a “perturbation” of the projected space, and cite smooth dependence of initial 
conditions; the trajectory going through a slightly perturbed initial point corresponding to 
such a case must enter V{0a+2) in a well-behaved way. In particular if zq belongs to Sq 
then a mapping from a neighborhood of zq on Sq to is well defined and Lipschitz- 

continuous. 

Suppose that case (ii) occurs. Then once an orbit of zq enters V (Oq,+i), it starts to manifest 
the dissipative behavior. In particular, if d{x, W^) < e, then after passing through V(Oa+i), 
d{x, < Ce'^, by Theorem 14.II Once the representative point leaves V{Oa+i)-, we may 

apply (i), viewing its position after leaving the neighborhood as an initial condition that 
does not re-enter the neighborhood V{Oa+i)- Thus when the representative point finally 
enters V{ 0 a+ 2 ), its distance to T has been contracted by an order of where 

1 < I'a+i — Ga+i and C absorbs both the constant C from Theorem 14.11 and a Lipschitz 
constant. 

Suppose now that zq G V(Ok), where k is now hxed. Then the trajectory leaves V{Ok) 
without increasing its distance from the unstable manifold, and passes into V{Ok+i) as just 
described. We may then apply Theorem 14.11 As the trajectory passes through V{Ok), its 
distance from the unstable manifold is contracted on an order of 
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Since the mapping contracts in the global dynamics and is Lipschitz (or contracting) in the 
local dynamics, simple inductions yields that as a representative point moves through the 
system, its distance from the manifold changes from e to to and 

so on. 

For a fixed i, the value r'j, representing a ratio of eigenvalues, is likewise fixed. The value 
is not; it depends on the distance between the representative point and the stable manifold 
as the trajectory enters 1^(0^), which changes from one instance to the next. For a given 
i, however, there is some maximal value that e* can take, since the system is constantly 
contracting towards the manifold and e* goes to 0 along with that distance. Thus there 
exists a global value, 1 < < i^i — Ci for all i and all e*, such that passing from the first to 

the p-th unstable manifold is a contraction of order . □ 


4.3 Existence of Parameter Sets 

Throughout the paper, we have put a number of restrictions on the parameters of the 
system. One must ask whether the specified inequalities may be satisfied. 

Lemma 4.2. There are sets of positive parameters values {ai}, i = l,...,n, and {pjk}, 
j,k = l,...,n, with non-empty interior such that the inequalities 0, © and (USD 
hold. 


Proof. First note that given {crj} the inequalities © and © are completely uncoupled and 
all trivially have positive solutions {pjk}- For each k and i = 1,2 these are: 

Pjk < ^ (20) 

CTk 

We note that with these inequalities the stable eigenvalues may be freely chosen to take 
any negative value and the unstable eigenvalues any positive values. Since the restrictions 
© concern only relative orderings of those eigenvalues at each Ok {k fixed), it is clear that 
© is satisfied for open subsets of the previously chosen sets. 

The final remaining inequality ([1© concerns on pk+i,k and Pk+ 2 ,k and so to finish the proof 
we need only to consider whether this restriction on those values is consistent with the 
previous restrictions on those parameters, namely © and ©, but not ©. Note that the 
inequalities in © can be written as: 

O'fc+l ^ , crA:+2 ^ 

- > (jfc and - > CTfe. 

Pk-\-l,k P/c+2,fc 

The constraint ()15p only requires that 

O'k+l ^ ^k +2 

Pk-\-l,k Pk-\- 2 ,k 

There is clearly no inconsistency in these inequalities. The final inequalities © involve 
each pk+i,k and pk+ 2 ,k independently of the others. First they require that for each k and 
f = l,2 

^k Pk+l,k^k ^ ^k- 
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This is equivalent to 




1 < Pk+l,k 


which can clearly be satisfied along with ()20p . If {p^+i^k}: i = 1, 2 have already been chosen 
to satisfy ([3]) and (fT^ . then values for pjk can be chosen to satisfy ([5]) by simply choosing 
them sufficiently large, i.e., for each j, k, j k,k + l,k + 2, pj^ must satisfy: 

^k+l Pk+i,k^k ^ Pjk^k j * — f)2. 


This can be rewritten as: 

^ c^k+i + • _ 1 n 

Pjk ^ Pk+i,kj * J-) 

Thus Pjk can take any value greater than the maximum of these two values. Previously we 
had only required (in (jl])) that these parameters satisfy: 


Pjk ^ 


(T,- 




Thus, with all other choices of consistent parameters, any large enough pjk, will also be 
consistent. □ 


For example if <7* = <7 are all the same, then ([3]), ([4]) and (fTKI) are satisfied if simple : 

0 < Pk+2,k < Pk+l,k < 1 


and 


Pjk > 1, for j k,k + l,k + 2. 

The dissipative requirement ([5]) will be satisfied if further: 

2(7 - pk+i,k < Pjk, for z = 1,2 and j k,k + l,k + 2. 


For instance, the parameter values: ai = 1, Pk+i,k = -9) Pk+ 2 ,k = -8 and pjk = 1.3 for 
j k,k + l,k + 2, (indices mod p) strictly satisfy all of the inequalities. 


5 Conclusion 


We proved in the paper that under some conditions the generalized Lotka-Volterra system 
admits a two-dimensional attractor that consists of saddles and the unstable manifolds 
joining them into a heteroclinic system. Thus, trajectories inside the attractor manifest 
completely regular features. However, behavior of wandering trajectories in the basin of the 
attractor could be treated as weakly chaotic. Indeed, one can introduce an oriented graph 
with vertices identified with the saddle equilibrium points Ok and edges identified with 
heteroclinic trajectories joining Ok and Ofc+i (belonging to the coordinate plane Pk,k+i) or 
Ok and Ok +2 (belonging to the plane Pk,k+ 2 )- It is possible to show that for each finite 
path through this graph there exists an open set of initial points in the basin such that 
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the trajectory going through any of these points follows the corresponding heteroclinic 
trajectories. The number of paths grows exponentially with the length, so the number of 
pieces of trajectories with different behavior (in fact all of them are (e, T) separated for 
some values of e and T) grows as T —>■ oo, the metric complexity function grows with time. 
A similar effect has been observed in [T], where it was called weak transient chaos. We 
intend to describe its properties in another publication. 

In the context of the motivating application, functional sequential dynamics in neural net¬ 
works, the two-dimensional heteroclinic attractor T in the phase space of dynamical system 
pi]) may be thought of as a mathematical image of diverse sequential dynamics based on 
the parallel performance of not one but two different modalities. Many interesting applica¬ 
tions of this may be found in cognitive science. For example, the learning and performing of 
sensory-motor human behaviors in many situations require the integration or binding of the 
sequential stimuli of one modality with the sequential stimuli of another. This seems to be 
the case in one of the most important cognitive functions: sequential working memory. In 
performance of a cross-modal working memory task, there may be two sequentially discrete 
neural processes (different chains of metastable states) that represent simultaneous neural 
activities corresponding to cross-modal transfer of information in the working memory (see 
e.g. [20]). 
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A The Topological Form of F 


Theorem ED states that F depends on the parity of p; this is because the number of 
connected components in the boundary of F depends on the parity of p. 

Proposition A.l. If p is even, then the boundary cIF has two connected components. If p 
is odd, then dT has one component. 


Proof. If p is even, then one connected component of the boundary will include the trajec¬ 
tories connecting the saddles On where n is even, and another will include the trajectories 
connecting the saddles Om, where m is odd. 

If p is odd, then as in the previous case, the saddles On where n is even are contained in 
the same component, and the saddles On where n is odd are likewise contained in a single 
component. Furthermore, Op-i and Oi are contained in the same component, where p — 1 
is even, since Oi = Op+i. Thus all saddles On are contained in the same component. □ 


Figure El will help visualize the situation. 
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Figure 5: A diagram representing p = 6 (left) and p = 5 (right). We clearly observe two 
distinct boundaries for p = 6, and one for p = 5. 


It appears from the diagram that when p is even, F is a cylinder, and when p is odd, F 
is a Mobius strip. We will formalize that intuition. We first observe that the components 
whose union is the unstable manifold F are not literally triangles, as they are depicted in 
Figure[5l but rather curved surfaces. We cite a definition from algebraic topology (e.g. [8]). 

Definition A.2. Suppose X is a compact Hausdorff space. A curved triangle in X is 
a subspace A of X and a homeomorphism H : T ^ A where T is a closed triangular 
region in the plane. A triangulation of X is a collection of curved triangles Ai, ..., in X 
whose union is X and such that for i ^ j, the intersection Ai D Aj is either empty, or a 
vertex of both Ai and Aj, or an edge of both. We also require, if hi is the homeomorphism 
associated with Ai, that when AjHAj = e is an edge of both, then the map hj^hi is a linear 
homeomorphism of the edge hf^{e) of Ti with the edge h~^{e) ofTj. 


All compact surfaces have a triangulation, but we will prove in particular that the decom¬ 
position of F into equilibria and the closure of their unstable manifolds, each restricted to 
the first orthant, forms a triangulation of F. 

Definition A.3. Consider the system of differential equations ©• For each ol, we call the 
closure ofWf^UOa a heteroclinic triangle, Ta. 


Note that: ^ 

Ta = W^UOaC Oa+l U 0^+2 U r«+i,„+2 

where rQ,_|_i^Q +2 is the heteroclinic orbit from Oa+i to Oo+ 2 - This follows from the invariance 
of the (a, a -|- 1, a -|- 2)-plane and Theorem 13.121 

Theorem A.4. A heteroclinic triangle is homeomorphic to a closed triangle in the 
plane. 


Certainly the boundary of a heteroclinic triangle T^, that is Oq, Oo+i, and O 0 + 2 , and 
the smooth paths connecting them, is homeomorphic to the boundary of a closed triangle 
in the plane. Further, the interior of is W^, by the unstable manifold theorem and 
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Theorem 13.121 homeomorphic to an open disk in the plane (and so also to the interior of 
a triangular region). The only complication may arise from the behavior of the unstable 
manifold near the edges or vertices. For instance, the could be “folded” as it approaches 
the edge connecting Oa+i and Oa +2 so that a neighborhood of a point on the edge is not 
locally homeomorphic to a point on an edge of a closed triangle in the plane. 

Consider the dynamics projected onto three consecutive dimensions; for simplicity of nota¬ 
tion, we will use the standard x, y, z coordinates, and label the saddle points on each axis 
Ox, Oy, and Oz, where there are heteroclinic connections Ox ^ Oy ^ Oz and Ox ^ Oz, 
which we denote as Txy, Vxz, and Tyz- We denote the corresponding heteroclinic triangle 
by Tx- 

Proposition A.5. Each orbit in W^{Ox) is uniquely identified with an angle 0 < < n. 


Proof. We define the usual ^-neighborhoods Ns{Ox), Ns{Oy), and NfiOz). By the Unstable 
Manifold Theorem, we may choose 5 small enough that W^{Ox) is the graph of a function 
of {y,z) in Ns{Ok). Using the monotonicity of the z coordinate inside the invariant region, 
we may also choose 5 small enough such that once a trajectory leaves N^{Ox) it cannot 
return to it. By the previous observations, the other coordinates (other than x, y, and z) 
of W'^{Ox) are all zero. 

Now consider that each orbit in W'^{Ox) has a unique intersection point w with the bound¬ 
ary of Ns{Ox). Since W'^{Ox) is a graph over {y,z) consider the projection w of w onto 
{y, 2 )-plane. Let cj) denote the angle of the ray through the origin and w with the z-axis. 
Note that 4> lies between 0 and 7 r/ 2 . The extreme angles 0 = 0 and 0 = 7 r /2 correspond to 
the heteroclinic orbits Txz and Txy respectively. □ 

Definition A.6. For a point p G W'^(Ox) U Ox define d{p) to he the distance from Ox to 
p along the orbit containing p. For a point p G W^{Oz) U Oz define e{p) to be the distance 
from p to Oz along the orbit containing p. 

Proposition A.7. The arc lengths d and e are finite and continuous where defined. 


Proof. By basic existence theory the orbits are smooth and therefore arc length is locally 
well-defined on them. 

Denote the solution with initial value p as p{t). Recall from the Unstable Manifold Theorem 
that p{t) —>■ Ox as f —)■ —oo and in fact 

\Ox - p{t)\ < for t < 0 , ( 21 ) 

where is the leading unstable eigenvalue, i.e. in this case the minimum of 

<72 —P 21 < 7 i and <73 —PsifJi, 

and e“ may be chosen to be arbitrarily small. 
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Now consider the length of the orbit: 


fO _ 

d{p) = / _|_ 2;2 

J —OO 

= f 1^x2 _|_ y2 _|_ 2-2 (It + f + 7/2 + df 

Jr J —OO 


( 22 ) 


If we substitute the equations dH) into the integrals, then substitute (f 2 T]) into the second 
integral one easily sees that this integral goes to zero as r —>■ —oo. Thus d{p) is hnite. 


Now consider the orbits starting at points q close to p. Fix e > 0 and chose r so that 
the remainder integral above is less than e/3. Since (f 2 TI) implies that |g(f) — p{t)\ goes to 
zero exponentially as t —>■ — oo, we can make the difference in the remainders less than 
2e/3 if (/(r) is sufficiently close to p(r). This we can accomplish by requiring q(0) and p(0) 
sufficiently close (by continuous dependence on initial conditions). We can also make the 
integrals from r to 0 less than e/3 by choosing g(0) close to p(0). For such q, \d{q)—d{p)\ < e. 

The arc length e{p) is finite for any p not in Ox U Txy U Oy since is a stable node in 
the xyz-subspace and all orbits approach Oz exponentially in time and so the same type 
estimates as above hold here. Continuity of this distance also follows by a similar proof as 
for d. □ 


Corollary A.8. All orbits on Tx have finite length and this length is a eontinuous funetion 
of an initial point on W'^{Ox) \ (ra,j^ U Oy U Fy^) 


This follows from the continuity of d and e where they are defined. 

Lemma A.9. The arc length d can be extended continuously to the edge Fy^. The arc 
length is finite and uniformly bounded for all orbits that make up Tx- 


Proof. Denote by the length of the orbit in W'^{Ox) identified by the angle 4> for 
0 < (f) < 7 r/ 2 . 

Note that distance is well defined along the edges r 2 ;y and Fy^. In fact since Fxy C W'^{Ox)0 
VF^(Oy), and Fy 2 C IF“(Oy) fl W^(Oz) the arguments above show that the lengths of these 
two solutions arcs are finite. Let Dxy denote the length of F^-y and Dyz the length of Fy^. 
Set 

d^T7l2 - hDxy T Dyz., 

i.e. the combined length of r 3 ;y and Fy^. 


For p S Fy 2 set: 

i.e. the length to p along r 2 ;y 

P G Fy^. 


d{p) = 0^/2 - e(/5), 

and Fy^. We claim that d thus dehned is continuous at 


Fix e > 0. Consider a <5 (sup norm) neighborhood N of Oy. For any sufficiently close to 
7 r/ 2 , the orbit on Tx indexed by passes through N. Denote by P~ the point where F^-y 
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intersects N and by P'*' the point where Ty^ intersects N. By continuity of d and e we may 
chose such a 6 sufficiently small that: 

— d{P~) < e/6 and Dyz — e{P^) < e/6. 

Further, require that 6 be sufficiently small so that: 5 'j < e/6. 

Given cj) sufficiently close to 7r/2, denote by Q~ and Q~^ the points where the orbit with 
angle cj) intersects N. Again by continuity of d and e for all (j) sufficiently close to n/2 we 
have: 

\d{Q-) - d{p-)\ < e/6 and e(Q+) - e(P+) < e/6. 


Let g be a point sufficiently close to p so that the above conditions on cj) are satisfied and 
such that \e{q) — e{p)\ < e/6. 

Now d{q) will be equal d{Q~) plus the arc length from Q~ to Q~^ plus the arc length from 
to q. Denote these arc lengths by and i{Q~^,q) respectively. 


Within the neighborhood N we can transform coordinates to straighten the unstable man¬ 
ifold Vyz- In this subspace this unstable manifold is the graph of a function of 2 ; only and 
has that has the form: 

y = ay + h{z). 

We thus straighten W‘^{Oy) by the coordinate change Y = y — Cy + h{z). 

By the Shilnikov variables technique m, we have that the stable variables satisfy: 

|(x(t),F(t))| < 0<t<r, 

while the unstable variable z satisfies: 

|z(f)| < 0<f<r, 


where r is the time of passage through N. We then obtain: 




< 

< 


f x/+ z^ dt 

/o 


-|- \z\ dt 


f \x\ + \y\ 

Jo 

Jo 


< 


6 


2 


+ 


1 

A“ + e“ 


< 


e 

6 ' 


(23) 
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We have that 


d{q) - d{p) = d{Q ) + 1{Q -.Q^) + .q) - {D^y + Dy^ - e{p)) 

= d{Q-) + ^(Q-, Q+) + e(g+) - e{q) 

- D,y + d{p-) - d{p-) - Dy, + e(P+) - e(P+) + e{p) (24) 
= d{Q-) - d{p-) + £iQ-,Q+) + e(g+) - e(P+) 

+ d{P~) - D^y + e(P+) - Dy^ + e{p) - e{q). 

Combining the above estimates we have: \d{q) — d{p)\ < e. 

Now let 0 <(/)<'7r/2 be the arc length of the orbit determined by the angle (f). By the 
continuity of d on all of Tx, except Ox and the continuity of e in a neighborhood of D^f, 
depends continuously on cj). It is thus uniformly bounded. □ 


Proof of Theorem \A.4\ For each coordinate pair p = {4>{p), d{p)), we normalize the ar- 
clength d by leting u = d/D^p Note that in this normalized distance u{Oz) = 1 and so u is 
continuous on Tx- 

We will show that the heteroclinic triangle Tx with corners , Oy, and Oz is homeomorphic 
to the triangle ABC in the plane with corners at A = (0,0), B = ( 6 , 1 / 2 ) and C = (1,0). 

Now for a given point w on T identified by (</>, u), consider the map: 

H ■. w ^ {u,v) = {u, h{u, 4>)). 
where v = h{u, 4>) is given by the continuous map: 

J ^ tan((/)/ 2 ) if 0 < u < 6 
1 2 li-b) tan((/)/ 2 ) ifb<u<l. 

The map if is a homeomorphism. It is one-to-one at all points. It maps Ox, Oy and Oz 
onto A, B and C respectively. At interior points it is a local homeomorphism because it 
is a composition of local homeomorphisms. At the corners Ox and Oz, considering the 
mappings as polar coordinates make it clear that if is a local homeomorphism there. It is 
also clear that H is one-to-one along the edges Vxz, ^xy and Tyz- □ 

Thus Figure [5] represents not merely an easy-to-understand representation of F, but a tri¬ 
angulation. We now investigate the orientibility of F. The most common definition of 
orientibility, in terms of normal vectors, is not useful in this situation, but having triangu¬ 
lated F, we may use instead a less common, but still standard definition (see e.g. m)- 

Definition A. 10. Consider some arbitrary triangle in the triangulation of a manifold. 
Assign to each triangle in the triangulation a value of clockwise. This assigns corresponding 
directions to each side of each of the triangles. Now let Aj and Aj be triangles sharing a 
side; observe that the common side has been assigned two different directions, i.e. two 
adjacent triangles with the same orientation conflict on their shared side. If this can be 
done without contradiction, then the surface is orientable. If a contradiction is reached, the 
surface is non-orientable. 
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It is a standard result of algebraic topology that orientibility is independent of the specific 
triangulation used. 

Using this definition and the specific triangulation we have defined, the following lemma 
can be easily proven. 

Lemma A.11. If p is odd, T is non-orientable. If p is even, T is orientahle. 


Proof. Let p be odd. We consider the triangulation by heteroclinic triangles. Start with 
Ai, the triangle defined by Oi, O 2 , and O 3 and without loss of generality assign it the 
clockwise orientation. Thus the direction on the Oi — O 2 side of the triangle is given the 
direction O 2 ^ Oi. Likewise, the triangle Aj defined by Op, Oi, and O 2 is assigned the 
clockwise orientation, and the Oi — O 2 side of the triangle is given the direction O 2 —>■ Oi. 
But Aj and Aj are adjacent triangles, and the fact that they give the same direction to the 
Oi — O 2 side means that the surface is non-orientable. 

Let p be even. We consider the triangulation by heteroclinic triangles. Giving the triangle 
dehned by Oi, O 2 , and O 3 the clockwise orientation, and thus giving the Oi — O 2 edge the 
O 2 Oi direction, we clearly observe that making the Op, Oi, O 2 triangle clockwise gives 
the Oi — O 2 boundary, the only place where orientibility could be broken, the conflicting 
Oi —> O 2 direction. □ 


We cite one more result, a classification theorem [ 8 ]. 

Theorem A.12. Given a compact connected triangulable 2-manifold Y with boundary such 
that dY has k components, Y is homeomorphic to X-with-k-holes, where X is S'^ or the 
n-fold torus Tn or the m-fold projective plane Pm- 


A “hole” in the sense of the theorem is a set homeomorphic to an small e open ball. 

Our extensive build-up makes the proof of the theorem almost trivial. 

Theorem A.13. If p is odd, then T is homeomorphic to a Mobius strip. If p is even, T is 
homeomorphic to a cylinder. 


Proof. Let p be odd. Of the possible homeomorphic images named in Theorem IA.121 only 
the projective plane is non-orientable; since we know from Lemma lA.ll that T has only 
one boundary component, T is homeomorphic to the projective plane with a ball removed, 
which is homeomorphic to the Mobius strip. 

Let p be even. Then T is a rectangle with two of its edges identified. Since the “right-hand” 
point of the bottom edge is identified with the “right-hand” point of the top edge, it is a 
cylinder. □ 
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